Using isosbestic points to extract interactions from structure factors 
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Inverting scattering experiments to obtain effective interparticle interactions for particles in solu- 
tion is generally a poorly conditioned problem. More accurate potentials can be obtained through 
the use of isosbestic points, values of k where the scattering intensity I{k) or the structure factor 
S{k) is invariant to changes in potential well-depth. These points also suggest a new extended cor- 
responding states principle for particles in solution based on the particle density or packing fraction, 
the second osmotic virial coefficient, and a new measure of potential range. 
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Extracting effective interparticle interactions from ex- 
periments is a central objective in chemical physics, be- 
cause these interactions govern physical behaviour [Q. 
But like many such inverse problems, this task is compli- 
cated. Experimental data is not perfect: statistical fluc- 
tuations occur, and results may not be obtainable over a 
complete parameter range. Moreover, the inversion pro- 
cedures are often ill conditioned: small differences in the 
initial input can result in large changes in the flnal out- 
put. This letter concentrates on a well-known class of 
such problems, namely the inference of (spherically sym- 
metric) effective interparticle pair potentials v'^^f (r) from 
experimental scattering data for particles in solution. For 
that reason lower effective volume fractions rj are inves- 
tigated than in the better studied case of simple liquids 
near the triple point, where the dominance of entropic 
hard-core forces]^ makes it difficult to distinguish be- 
tween different attractive potentials!^, ^. Even though 
at lower volume fractions attractive interactions can de- 
termine the structure^, this letter shows that inversions 
are equally problematic, albeit for different reasons. 

For particles in solution, the scattering intensity I{k), 
measured by light. X-rays, or neutrons, is usually in- 
terpreted by dividing I{k) by the single particle form 
factor P{k), to obtain the structure factor S{k) = 
I{k)/{pP{k))^; here p = N/V is the density. From 
S{k), inferences are then made about the form of the 
inter-particle interactions v'^^^ {r). However, this inver- 
sion is ill-conditioned: Within typical experimental ac- 
curacy, a single S{k) at a given state-point can be inter- 
preted with a wide range of different v''^-^{r)'s. But, as 
this letter demonstrates, these problems can be circum- 
vented by introducing a new concept, isosbestic points 
- values of the wave-vector k for which S{k) is invariant 
under changes of the attractive potential well depth. Isos- 
bestic points determine the effective range of the v'^ffir), 
and also suggest an extended corresponding states prin- 
ciple for particles in solution, described by three experi- 
mentally accessible variables: the range defined by these 
points, the particle density p and the reduced second os- 
motic virial coefficient = ^2/^2^'^, where is the 
second virial coefficient of a hard-sphere (HS) system. 
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FIG. 1: A single set of experimental data can be interpreted 
by many different effective potentials, (a) The experimental 
SANS data (circles) on a microemulsion at packing fraction 
rj = 0.075 is from [6|. The open triangles denote the S{k) from 
the Baxter model] 7 1 for — —1.60, while the other S{k) are 
calculated for the potentials shown in (b). These include SWs 
with Asw = 0.02 (solid lines) and Asw ~ 0.2 (dashed), an 
AO form with Aao = 0.4 (long-dashed), a generalised deple- 
tion form with A_d = 0.2(dotted) and a LJ-12 (dot-dashed) 
potential. The line-styles in the two graphs correspond with 
each other. 

Fig. |l| helps set the stage and introduce the problem 
to be addressed here. Experimental data]^, taken from 
small-angle neutron-scattering (SANS) of a microemul- 
sion composed of small water droplets, coated with a 
layer of surfactant (AOT), is compared to some theo- 
retical structure factors S{k). In the original letter|^ 
a best fit, using an approximate integral equation tech- 
nique, yielded a hard-core diameter of cr = 60 A, a pack- 
ing fraction 77 = na^/d — 0.075, and an attractive square 
well (SW) potential with a range Asw = 0.02cr and a 
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well depth of (3v{a) = 3.85. The quality of the inte- 
gral equation fit was then confirmed by independent com- 
puter simulations. However, as strikingly demonstrated 
by Fig. ^ a wide range of other potentials, depicted in 
Fig. |l|(b), also lead to fits of similar accuracy. These 
include HSs with attractive SWs of two different ranges 
Asw, an Asakura-Oosawa (AO) depletion potential[|| of 
range Aao i an alternate simplified depletion potential of 
range showing the effects of solvation layers]^, as 
well as a Lennard-Jones-n (LJ-n) potential, defined as: 

•M-<(7)""-(7)")' 

with n = 12. The S{k) are calculated within the Percus 
Yevick (PY) integral equation closure |l^ , and are close 
to the analytical PY solution of Baxter's modelQ. 

Improving significantly on these pioneering experi- 
ments is hard. iS'(fc) is obtained by dividing I(k) through 
by P{k), which is only approximately determined, and 
which rapidly decays to zero for increasing A:cr, so that 
achieving better accuracies for a larger ka range is very 
difficult. Moreover, there are few alternative methods to 
directly measure the interparticle interactions. In con- 
trast, for simpler atomic or molecular systems, gas-phase 
experiments or ab-initio calculations can provide accu- 
rate independent estimates of the interactions, and im- 
pose important constraints on the potential forms with 
which to fit the experimental data||ll|. Although in- 
creasingly accurate techniques being developed for the 
real-space measurement of the v'^^^{r) of particles in 
solution p^, these are often limited to certain particle 
types or sizes. 

The upshot of Fig. ^ is that deducing an effective po- 
tential from experimental scattering data at one state 
point is difficult; the inversion is ill-conditioned. Clearly 
more information is needed to interpret the datap3[. 

But first, what can be inferred from a single S{k)l One 
hint comes the Baxter model, which is completely deter- 
mined by the packing fraction 77 and the reduced osmotic 
virial coefficient i?2 or equivalently, the Baxter param- 
eter, defined as r = - B^). The S{k) in Fig. 0, 
with effective t parameters ranging from r — 0.093 to 
T = 0.10 are all close to the Baxter model result with 
T = 0.096. To first order therefore, measuring a single 
S{k) in this regime of low packing fraction, commonly 
encountered for particles in solution, does not result in 
much more information about the effective potential than 
the reduced second virial coefficient i?2- Just as in sim- 
ple liquids inversions are hard because the S{k) resemble 
those of HSs||^, here difficulties arise because they are 
close to the Baxter S{k). 

To make further progress, measurements at other state 
points are necessary. This is done in Figs, ^(a)-(d), which 
depict the S{k) calculated for four LJ-n potentials at dif- 
ferent temperatures. A best fit to the Baxter model S{k) 
at each temperature would result in B2 as a function of 



temperature. But there is clearly more information in 
these curves: Each of the four v'^^^{r) results in a differ- 
ent set of isosbestic points, where the S{k) is invariant 
for different temperatures. Fig. ||(e) shows that the first 
isosbestic point kia increases with increasing n, reflect- 
ing the decrease of the range of the LJ-n potential (|l|) 
depicted in Fig. ||(f). (Error bars in Fig |^(e) reflect the 
approximate nature of the isosbestic points.) 

These observations can be quite easily rationalised 
by the following simple theory: If the total correlation 
function h(r) = g{r) — 1 is split into two parts, with 
ho{r) = —1 for r < a, and hi{r) = g{r) — 1 for r > cr, 
then the structure factor S{k) = 1 + ph{k) simplifies to 

s{k) = i + ^n{k) + pU{k), (2) 

where ji(fc) is the first spherical Bessel function, and 
hi{k) is the Fourier Transform (FT) of hi{r). For po- 
tentials with a hard-core, such as most of those depicted 
in Fig. |l|(b), this is exact, but even for the LJ-n potentials 
this is a good approximation |Q. Since /lo(fc) = ^f^ji(fc) 
is independent of temperature (barring small effective a 
effects), Eq. | suggests that isosbestic points occur when- 
ever hi{k) = 0. 

It has already been pointed out that in many exper- 
imentally relevant cases with low rj, the gir) are sur- 
prisingly well approximated by a simple form g[r) = 
exp(— /3w(r))|^, [l5|, as explicitly demonstrated in Fig. ||. 
This simple approximation, equivalent to taking a Mayer 
function|^] for h{r), works best for particles with short 
range attractive potentials. These can become quite 
deep, leading to large values of g{r) near contact, well 
before the system crosses a liquid-liquid or liquid-solid 
phase-line. Within this approximation, the dominant ef- 
fect of varying the temperature is to change the ampli- 
tude of hi{r) (as demonstrated in Fig. §(b)). To first 
order, the period of hi{k) is not affected and each ka 
where hi{k) — leads to an isosbestic point. For an 
infinitely narrow potential, the isosbestic points would 
be at fcjifj = riTT, but for a finite potential range there 
is a phase-shift, which explains why the kncr move pro- 
gressively further from nir with increasing range. The 
excellent accuracy of this simple Mayer function theory 
for the first isosbestic point kia is demonstrated by the 
solid line in Fig. ||(e). 

The examples depicted in Fig. |^ are at a relatively 
low packing fraction, but the isosbestic points are robust 
up to packing fractions of at least 77 = 0.2, as shown in 
Fig. ^(a)[|6|. These points are accurately described by 
the simple theory, described above. 

If varying a parameter leads to isosbestic points, then 
this suggests that the well-depth of v*^^^ {r) is changing 
while the range is not. When combined with the varia- 
tion of i?2, this can help fix the form of v'^^^ {r). How- 
ever, this information may still not be sufficient, since 
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FIG. 2: (a)-(d) S(k) were calculated at 77 = 0.0576 and at 
different temperatures for each of the four LJ-n potentials 
v^^ ^ (r) given by Eq. (^). Different v''^^{r) lead to different 
isosbestic points fe„o-. (e) The first isosbestic point kia ap- 
proaches 77 with increasing n. A simple theory (solid line), 
described in the text, fits the data well, (f) The range of the 
LJ-n potentials in (a)-(d) decreases with increasing n. 
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FIG. 3: (a) The PY approximation (solid line) very accu- 
rately reproduces these molecular dynamics simulations p^] 
of a LJ-12 potential at /3e = 1.6, and rj = 0.1. The simpler 
g{r) — exp{—j3v{r)) form (dashed lines) is also accurate, (b) 
When the temperature is changed, the dominant effect on 
g{r), calculated here with PY, is to change its amplitude. 

the well-depth of v'^^^ [r) may depend in an (unknown) 
non-linear fashion on the parameters like the tempera- 
ture T, the pH, and salt or other additive concentration. 
Furthermore, as shown in Fig. |^(b), different potential 
shapes, picked to have the same isosbestic points, can 
generate similar S{k) for each value of 

This similarity of the >S'(fc) ties in with the work of 
Noro and Frenkel (NF), who proposed that many prop- 
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FIG. 4: Fig (a) The position of isosbestic points barely varies 
with density for LJ-12, shown here for three temperatures 
/3e = 1.5,1.33,1.09 at each density r) = 0.05,0.1,0.15,0.2. 
The line with the arrow denotes feicr for increasing 77. (b) 
Three potentials leading to the same isosbestic points: SW 
(Asw = 0.26) (solid lines), AO (Aao = 0.6) (dashed 
lines) and LJ-22 (dot-dashed lines) for r) — 0.05 and B2 = 
0, —0.719, —1 (plots in descending order of B2 a,t k — 0). 

erties of particles in solution could be understood from an 
an extended corresponding state principle based on the 
variables 77, B2, and an effective well depth /3eQ. How- 
ever, the three potentials shown in Fig. ^(b) do not have 
the same well-depth, while still showing similar S{k). 
This suggests an alternative corresponding states prin- 
ciple based on the variables rj, B2, and an effective range 
Ae// related to the isosbestic points. In fact, as NF 
point out, the potential range is not always uniquely 
defined when comparing different w'^-^'^(r)'s. They sug- 
gested a non-linear mapping based on B2 to derive an 
effective SW range for each potential they consider. I 
used a simpler linear mapping to define the range A^ff 
of a given potential as identical to that of a SW with 
the same fcicr. For the LJ-n potentials this criterion 
reduces to the distance r — a where Pv{r) reaches 1/4 
of its maximum depth, for the AO potential the effec- 
tive SW range is Ag// ~ 0.42A^o, and for the hard- 
core Yukawa potential, with an attractive tail of the 
form v{r) = eexp(— r/A) for r > a, the mapping is 
Ae// ~ 1.15A. The results are depicted in Fig. ||, and 
show that these simple linear mappings successfully de- 
fine an effective range for each potential type. An accu- 
rate approximation of kia for the (effective) SW is given 
by kia w 7r/(l + Ae///2 -f Agyj,/12), which follows from 
an expansion of the exact solution to the simple theory 
of isosbestic points. 

Another argument for including the ^eff (instead of 
the well-depth) as one of the relevant 3 variables is 
that the range determines the topology of the phase- 
diagram jl^. The critical range below which the fluid- 
fluid transition becomes metastable to the fluid-solid 
line is at Ag// « 0.15 for all four potentials (SW, AO, 
Yukawa and LJ-n), something consistent with what NF 
found. In addition, A^jf is directly experimentally ac- 
cessible through kia. For example, taking advantage of 
the proximity of a metastable fluid-fluid critical point to 
enhance nucleation, a proposal recently made for pro- 
tein solutions lis], would involve tailoring solution con- 
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FIG. 5: Defining effective range parameters from isosbestic 
points. For each potential type a simple linear mapping to an 
effective square well range Ae/ / was used (as described in the 
text). The long-dashed line denotes an accurate analytical 
approximation for the isosbestic points of square-well fluids 
which can be used to estimate the Ae/ / for different potentials 
directly from the kia. 

ditions such that the measured isosbestic point is just 
above fcicr = 2.92, so that Ae// is just below 0.15. 

In summary then, extracting effective potentials 
y'^ff{r) from S{k) at the lower packing fractions typically 
encountered for particles in solution is a poorly condi- 
tioned problem. Whereas in simple liquids the HS model 
describes the dominant features of the S'(fc), here the 
Baxter model appears to be the fundamental underlying 
model for S{k) around which different attractive poten- 
tials only induce a mild perturbation. More information 
can be obtained by studying the isosbestic points, which 
help determine the range Ae// of the potentials. To- 
gether with rj and B2 , Ae/ / can be used to define an ex- 
tended corresponding states principle. Particles in solu- 
tion with similar values of these three parameters should 
have similar properties, such as the relative stability of 
the fluid-fluid and fluid-solid binodals. In other words, 
many solution properties could be deduced without hav- 
ing to actually invert to an explicit form of v^-^^{r). 

Another advantage of a description based on these vari- 
ables is that they are directly experimentally accessible. 
In fact, numerous examples of isosbestic points can be 
found in the literature, ranging from colloids with short- 
range sticky coats pjt, to colloid-polymer mixtures [pO|, 
to globular proteins to magnetic colloids |Q and to 
colloid-micelle mixtures [p3|. Interestingly, these points 
should also appear in I{k), although determining Ae// 
still depends on an accurate P{k). Future work will con- 
sider the role of size polydispersity (expected to mainly 
affect higher order A:„cr) and the existence of isosbestic 
points for non-spherical objects. 

The author thanks P. Bartlett for first bringing isos- 
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